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Abstract 

We explicitly construct global strict Lyapunov functions for rapidly time-varying nonlinear control systems. The Lyapunov 
functions we construct are expressed in terms of oftentimes more readily available Lyapunov functions for the limiting dynamics 
which we assume are uniformly globally asymptotically stable. This leads to new sufficient conditions for uniform global 
exponential, uniform global asymptotic, and input-to-state stability of fast time-varying dynamics. We also construct strict 
Lyapunov functions for our systems using a strictification approach. We illustrate our results using several examples. 
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1 Introduction 

The stabilization of nonlinear and nonautonomous con- 
■ trol systems, and the construction of their Lyapunov 
functions, are challenging problems that are of signifi- 
cant ongoing interest in the c ontext of robustness analy- 
sis a nd contro ller design; see jMalisoff fc Mazend i|2005j) 
and iMazenc fc Bowong (2004). One popular approach 
to guaranteeing stability of nonautonomous systems is 
the averaging method in which exponential stability of 
an appropriate autonomous system implies exponential 
stability of the original dynamics wh en its time variation 
is sufficiently fast. See Khalil (2002) for related results. 

i 

The preceding results were extended to more general 
rapidly time- varying systems of the form 

x = f(x,t,at), x e K™, t e K, a > (1) 

in lPeuteman fc Aevelsl l)2002f) . where uniform (local) ex- 
ponential stability of (1) was proven for large values of 
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the constant a > 0, assuming a suitable limiting dynam- 
ics 

x = f(x,t) (2) 

for (1) is uniformly exponentially stable. (We specify the 
choice of / in our main theore m and examp les below.) 
This generalized a result from (Hal£ 119801 pp. 190-5) 
on a class of systems (1) satisfying certain periodicity or 
almost periodicity conditio ns. The main arguments of 
iPeuteman fc Aevelsl ()2002l) use (partial) averaging but 
do not lead to explicit Lyapunov functions for (1). 

In this work, we pursue a very different approach. In- 
stead of averaging, we explicitly construct a family of 
Lyapunov functions for (1) in terms of more readily avail- 
able Lyapunov functions for (2), w hich we again assume 
is asy mptotically stable. While IPeuteman fc Aevelsl 
(2002) assumes (2) is uniformly locally exponentially 
stable, we consider the case where (2) is assumed uni- 
formly globally asymptotically stable (UGAS), in which 
case our conclusion is that (1) is UGAS (but not nec- 
essarily exponentially stable) when a > is sufficiently 
large. While global exponential and global asymptotic 
stabilities are equivalent for autonomous systems un- 
der a coordinate change in certain dimensions, the 
coordinate changes are not explicit and so do not lend 
the mselves to explicit L yapunov function constructions; 
see lGrune et all p999|) . The importance of the problem 
we consider lies in the ubiquity of rapidly time-varying 
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systems in practical applications (e.g., suspended pen- 
dulums subject to vertical vibrations of small amplitude 
and h igh frequency , Raleigh's and Duffing's equations 
from (Khalil, 2002 1 Chatper 10), and s ystem s arising 
in identification in beuteman fe Aevelsl |2002|) and be- 
low) and the essentialness of Lyapunov functions in 
robustness analysis and controller design for these ap- 
plications. 

In particular, we show that ass umptions similar to those 
of (jPeuteman fc Aevelsl l2002l Theorem 3) imply that 
(1) is uniformly globally exponentially stable; our Lya- 
punov function constructions are new even in this par- 
ti cular case and our results a re complementary to those 
of iPeuteman fe Aevelsl 1)20021) . The Lyapunov functions 
we construct are also input-to-state stable (ISS) or inte- 
gral ISS Lyapunov functions for 

x = f(x,t,at)+g(x,t,at)u (3) 

under appropriate conditions on / and g; see Remark 10. 

In Section 2, we provide the relevant definitions and 
lemmas. In Section 3.1, we present our main sufficient 
conditions for uniform global asymptotic and exponen- 
tial stability of (1), and for the stability of (3), in terms 
of limiting dynamics (2) . This theorem leads to explicit 
constructions of Lyapunov functions for (1) and (3) in 
terms of Lyapunov functions for (2), in Theorem 8. In 
Section 3.2, we provide an alternative Lyapunov func- 
tion construction theorem for (1) not involving any lim- 
iting dynamics. We prove our main results in Sections 
4 and 9. We illustrate our theorems in Sections 6 and 7 
using a friction control model and other examples. We 
close in Section 8 by summarizing our findings. 

2 Assumptions, definitions, and lemmas 

We study (1) (which includes dynamics (2) with no a de- 
pendence, as special cases) in which we always assume / 
is continuous in time feR:=(-oo, +oo), continuously 
differentiable (C 1 ) in a; € R", null at x = meaning 

f(0,t,at) = f(0,t) = VteK, a>0 (4) 

and forward complete, i.e., for each a > 0, x £ R™, 
and t £ R>o := [0, oo) there is a unique trajectory 
[to,oo) 3 t i ^ <fi(t;t ,x ) for (1) (depending in general 
on the constant a > 0) that satisfies x(t a ) = x a . We set 
N = {1, 2, 3, . . .} and let Z denote the set of all integers. 
We say N : R> — ► R>o is of class M and write N £ M 
provided 

lim r)N(r)) = 0. (5) 

77— >+oo 

A continuous function 6 : R>o — > R>o is positive def- 
inite provided it is zero only at zero. A positive defi- 
nite function 6 is of class IC (written 5 £ IC) provided 
it is strictly increasing; if in addition S is unbounded, 



then we say that 5 is of class /Coo and write <5 £ /Coo- A 
continuous function f3 : R>o X R>o — * R>o is of class 
ICC (written [3 £ ICC) provided (a) (3(-,t) £ /Coo for all 
t > 0, (b) [3(s, ■) is nonincreasing for all s > 0, and (c) 
for each s > 0, f3(s, t) — ► as t — ► +00. A positive def- 
inite function S is called o(s) provided 8(s)/s — > as 
s — > +00. We always assume there exists p £ /Coo such 
that \f(x,t, at)\ < p(\x\) everywhere. 

We next define our stability properties for (2). The same 
definitions apply for (1) for any choice of the constant 
a > 0. We call (2) uniformly globally asymptotically sta- 
ble (VGAS) provided there exists /3 £ ICC such that 

\(f>{t;t ,x )\ < [3(\x \,t — t a ) Vt>t > 0, x a £ K" (6) 

where | • | is the usual Euclidean norm and <p is the flow 
map for (2). We call (2) uniformly globally exponentially 
stable (UGES) provided there exist constants D > 1 and 
A > such that (6) is satisfied with the choice 

(3{s,t) = Dse- Xt . (7) 

The converse Lyapunov function theorem implies (2) is 
UGAS if and only if it has a (strict) Lyapunov function, 
i.e., a C 1 V : 1" x R> -> R> that admits 6%, 62 £ /Coo 
and S 3 £ IC such that for all t £ K> and £ £ K™, 
we have both (LI) ^fl^) < V(£,i) < <5 2 (|£|) and (L2) 
Vt(£,t) + V£({,t) ffat) < -foflgn, where the subscripts 
on V indicate partial gradients; see iBacciotti fc Rosier 
(2005). When (2) is UGES, the proof of JKhaiIT T2002. 
Theorem 4.14) shows: 

Lemma 1 Assume (2) satisfies the UGES condition 
(6)-(7) for some constants D > \_and A > and that 
there exists K > A such that \ (df/d£) (£, t) \ < K for 
all £ £ 1" and t £ M>o. Then (2) admits a Lyapunov 
function V and constants ci, C2, C3 > such that 

ci|£| 2 < V&t) < c 2 \e, \V&,t)\ < 0s|£|, 
and V t (Z,t) + Vt(U)f(t,t) < -l^ 2 1 ' 

hold for all t £ R> and £ £ R" . 

Definition 2 Given S £ IC, the dynamics (2) is said to 
be 5-compatible provided it admits a Lyapunov function 
V £ C 1 and two constants c £ (0, 1), c > such that: 

«Pa V t (S,t) + V&,t)f(£,t) < -cS 2 m VC.i. 

P2 W€,t)\ < 5(\$\) and \f(£,t)\ < 5(\Z\/2) V$,i. 

P 3 S(s) < cs Vs > 0. 

Remark 3 Note the asymmetry in the bounds on |V S | 
and l/l m P2. If (2) satisfies the assumptions of Lemma 
1, then it is ^-compatible with 5(s) = (03 + 2K)s. How- 
ever, by varying 8 (including cases where S is bounded), 
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one also finds a rich class of non-UGES (5-compatible dy- 
namics; see e.g. Section 6.1 below. 

We also consider the nonautonomous control system 

x = F(x,t,u) (9) 

which we always assume is continuous in all variables and 
C 1 in x with F(0, t, 0) = 0, and whose solution for a given 
control function u G U(:~al\ measurable locally essen- 
tially bounded functions [0, oo) — > R m ) and given initial 
condition x(t ) — x we denote by t i— * <^>(t; t , x , u). 
We always assume (9) is forward complete, i.e., all tra- 
jectories (/>(•; t , x a , u) so defined have domain [i , +oo). 
We next recall the input-to-state stable (ISS) a nd inte- 
gral in put- to-state stable (iISS) properties from lSontad 
(1989) and lSontael l)1998j) . Let |u|j denote the essential 
supremum of u G U restricted to any interval I C R> . 

Definition 4 (a) We say that (9) is ISS provided there 
exist 7 G /Coo o,nd f3 £ K,C for which 

\(j>{t;t ,x ,u)\ < f3{\x \,t-t )+j(\u\ [tott+to] ) (10) 

holds whent > t Q > 0, x £ R", cmdu G 14. If in addition 
[3 has the form (7), then we say that (9) is input-to-state 
exponentially stable (ISES). (b) We say that (9) is iISS 
provided there exist fi, 7 G /Coo and (3 G ICC such that 



Since (9) has an ISS Lyapunov funct i on wh en it is ISS 
(b y the arguments oflSontag fc Wanel 1)19951) 1. the proof 
of l)Angeli et all. I2000L Theorem 11 shows that if (9) is 
ISS and autonomous, then it is also iISS, but not con- 
versely, since e.g. x = — arctan(ir) + u is iISS but not 
ISS. The nex t lemma follows from th e argu ments used in 
lAneeli et al] l)2000D . lEdwards et al\ i|2000D . and lSontad 



Lemma 6 If (9) admits an ISS (resp., iISS) Lyapunov 
function, then it is ISS (resp., iISS). 

3 Statements and discussions of main results 

3.1 Main theorem and Lyapunov function construction 



We sh ow that hypotheses similar to those of ijPeuteman fc Aevelsl 
120021 Theorem 3) ensure that (1) is UGES. In fact, we 
show our conditions imply 



f(x, t, at) + u, x £ 



u £ 



(12) 



is ISS when a > is sufficiently large; see Remark 10 be- 
low for results on the more general systems (3). Our main 
assumption will be: There exist 6 £ K,, a (5-compatible 
dynamics (2), and N £ M. (cf. (5) above) such that for 
all x £ R™, all r £ R and sufficiently large n > 0, 



rt +t 

fi{\(j){t;to,Xo,u)\) < f3(\x \,t-t ) + J 7(|u(s)|)ds 



holds when t > t Q > 0, x £ 



and u G IA. 



A function V : R" x R>o — >• R>o is called uniformly pos- 
itive definite provided s 1— > vaf{V(x,t) : t > 0, \x\ = s} 
is positive definite in which case we write V £ UPD. We 
call V £ UPD uniformly proper and positive definite pro- 
vided there exist 61,62 £ /Coo such that condition (LI) 
above is satisfied in which case we write V £ UPPD. 
The following Lyapunov function notions agree with the 
usual ISS and iISS Lyapunov function definitions when 
(9) is autonomous, because functions x *= ^00 are in- 
vertible. Note that v in (11) need not be of class /C. 

Definition 5 Let V £ C^nUPPD. (a) We callV an ISS 
Lyapunov function for (9) provided there exist %, £3 G 
/Coo such that for all t £ R> , £ £ R n , and u £ R m : 
[\u\ < x(\Z\) => V t (S,t) + V&,t)F(S,t,u) < -6 3 (m- 
(b ) We call V an iISS Lyapunov function for (9) provided 
there exist A G /Coo om,d a positive definite function v : 
R>o — * R>o such that 

Vt(S,t) + V&,t)F(S,t,u) < -v(\£\) + A(\u\) (11) 

holds for all t £ R> 0; £ G R" , and u £ R m . 



{f{x,l,r l 2 l)-f{x,l)}dl 



<6{\x\/2)N{n) (13) 



which is similar to ijPeuteman fc Aevelsll2002l Property 
2) in the special case where 6(s) = 2s. Two more ad- 
vantages of our result are (a) it applies to cases where 
(2) is UGAS but not necessarily UGES (cf. Section 6.1 
below) and (b) its proof leads to explicit Lyapunov func- 
tions for (12) (cf. Theorem 8 below). See also Theorem 
9 below for cases where df/dx is not necessarily glob- 
ally bounded. Recall the definition of compatibility (in 
Definition 2) and the requirement N G M. in (5). 

Theorem 7 Consider a system (1). Assume there exist 
5 £ /C, a 8-compatible UGAS system (2), two constants 
rj > and K > 1, and N £ M such that (13) holds 
whenever n > rj , x £ R" and r £ R and such that: 



(x,t) 



< K 



%{x,t,at) 



\f{x,t,at)\ < 6(\x\/2) Vi 



< K 

[,x£l 



and 

, a > 0. 



(14) 



Then there is a constant a > such that for all constants 
a > a, (1) is UGAS and (12) is iISS. If in addition 
5 £ /Coo, then (12) is ISS for all constants a > a. In 
the special case where (2) is UGES, (1) is UGES for all 
constants a > a and (12) is ISES for all constants a > a. 
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By (4), the condition \ f(x,t,at)\ < S(\x\/2) in (14) is 
redundant when 5 has the form S(s) = fs for a constant 
f > 0, since f can always be enlarged. Our proof of 
Theorem 7 in Section 4 below will also show: 

Theorem 8 Let the hypotheses of Theorem 7 hold for 
some 5 £ JC, and V £ C 1 be a Lyapunov function for (2) 
satisfying the requirements of Definition 2. Then there 
exists a constant a > such that for all constants a > a, 



v U-^f SL-fe /.*{/(& J. «0 ^ /(£> 0} d/ da, t 



is a Lyapunov function for (1) and an USS Lyapunov 
function for (12). If also 5 £ /Coo, then V la] is also an 
ISS Lyapunov function for ( 12) for all constants a > a. 

3.2 Alternative result 

The proof of Theorem 7 constructs strict Lyapunov func- 
tions for (1) in terms of strict Lyapunov functions for the 
limiting dynamics (2). It is natural to inquire whether 
one can instead construct strict Lyapunov functions for 
(1) by strictifying nons t rict L yapunov functions for (1); 
see IMalisoff fc Mazencl l|2005f) where the strictification 
approach was applied to nonautonomous systems that 
are not rapidly time-varying. In this section, we extend 
this approach to cover (1). Our strictification result has 
the advantages in certain situations that (a) it does not 
require any knowledge of limiting dynamics, (b) it allows 
the derivative of the nonstrict Lyapunov function to be 
zero or even positive at some points, and (c) it does not 
require (1) to be globally Lipschitz in the state. In the 
rest of the section, we focus on systems with no controls, 
but the extension to the control system (3) for appropri- 
ate g can be done by similar arguments. We assume: 

H. There exist V £ C 1 n UPPD, W £ UPD, a C 1 
function 9 : M. n x R> — * K, a bounded continuous 
function p : K — > K, and constants c, T > such 
that for all x £ R", t > 0, a > 0, and k £ Z, we have: 

HI. V t + V x f(x, t, at) < -W(x, t) + p(cd)Q(x, t) 



r (fe+i)T, 

kT 1 

H3. V > c\Q\, W > cmax{|0|, \Q t +®xf(x, t, at)\} 



H2. p(r)dr = 



Here and in the sequel, we omit the argument (x, t) of 
V, W, O, and the partial gradients of V and O whenever 
this would not lead to confusion. In Section 5, we prove: 

Theorem 9 If Assumption H. holds, then there exists a 
constant a > such that for all constants a > a, 



UW(x,t) 



V{x,t) - (f/_ 1 (f s p(al)dl) ds) @(x,t) 



(15) 



is a Lyapunov function for (1). In particular, (1) is 
UGAS for all constants a > a. 

4 Proof of Theorems 7 and 8 and remarks 

To make our arguments easy to follow, we first outline 
our method for proving these theorems. First, we give 
the proof of Theorem 7 for the special case where (2) is 
UGAS and 6 £ /Coo which includes the proof of Theorem 
8 for the S £ /Coo case. Then we indicate the changes 
required if S 6 K is bounded. Finally, we specialize to 
the special case where S(s) — fs for some constant f > 
which will prove the ISES assertion of Theorem 7. 

We first assume that (2) is UGAS and 5 £ /Coo, and 
we prove the ISS property for (12) for large constants 
a > 0. In what follows, we assume all inequalities and 
equalities hold globally, unless otherwise indicated. Let 
r\ be as in the statement of the theorem, and fix a = ?y 2 
with n > rj , u £ U, and a trajectory x(t) for (12) and 
u, with arbitrary initial condition. Set 



z{t) = x{t) + R a (x(t),t), 



(16) 



where 



R a {x,t) = -- 



t-2/ri J s 



{f(x,l, V 2 l)- f(x,l)}dlds. 



This is well-defined since we are assuming our dynamics 
are forward complete. Set 

P (t,l) = f(x(t),l, V 2 l)-f(x(t),l). 

With p so defined, one ea sily checks (as was done e.g. in 
IMalisoff fc Mazend 1)20051) 1 that for any r > and t £ E, 



dtJt- 



p(t,l) dlds 



rp(t, t) ~ Ii r P{% I) dl + J J ^ (t, I) dl ds 



(17) 



and 



flrlsPfrl) dlds 



< 



t-T<Kt 



\p(t,l)\. (18) 



Taking r = 2/n, (17) multiplied through by —T)/2, 
f{x{t),t,r, 2 t)-p{t,t) = f{x{t),t), and (16) give 

z(t) = f{z{t),t) + (f(x(t),t) - f(z(t),t)) + u(t) 

+1 lU p(*> o dl - 1 \ ft— — ft (§£(*(*)> h v 2 i) (19) 

r) V rj \ 

-$£(x{t),lj) dlds} (f(x(t),t,r,H) + u(t)) . 
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Let V, Si, and S 2 satisfy the requirements P1-P3 from 
our compatibility condition (in Definition 2) and (LI). 
By Pi with £ = z(i) and (19), the derivative of V(z, t) 
along the time- varying map z(t) defined in (16) (which 
we denote simply by V in the sequel) satisfies 

V < -c6 2 (\z(t)\) + V s (z(t),t)(f(x(t),t)-f(z(t),t)) 

+%Vt(z(t),t)fl 2/v p(t,l)dl 

-W(z{t),t) 

x (f(x(t),t, i?t) + u(t)) + V t {z{t),t) u(t). 
We deduce from (13), (14), (16), and P 2 that 

V < -c5 2 (\z(t)\) + KS(\z(t)\)\x(t)-z(t)\ 

+^S(\z(t)\)\fl 2/ri p(t,l)dl\+S(\z(t)\)Ht)\ 
+ § (\f(x(t),t, V H)\ + |u(i)|)J(|s(t)|) 

*SL,J>k&i&s 

< -c8 2 (\z(t)\) + K6(\z(t)\)\R a (x(t),t)\ 

+%5(\z(t)\)N(r 1 )5(\x(t)\/2)+5(\z(t)\)\u(t)\ 

+$K5(\z(t)\){8(\x(t)\/2) + \u(t)\}- 



Moreover, (14), the definition of p, and P 2 give 

\R a (x(t),t)\ < |//_ 2/ J>(M)|dZd s 
< §<KW*)I/ 2 )- 



(20) 



Combining these inequalities and grouping terms gives 

V <-c6 2 (\z(t)\) + 6(\z(t)\)\u(t)\ 

+6(\z(t)\) (±K + %N{r,)) {S(\x(t)\/2) + |u(i)|} . 

On the other hand, (16), (20), and P 3 give 

W)\ > \x(t)\ - l\x(t)\ > ±\x(t)\ (21) 
when 77 > max{2c, r] }. Since 6 G JC, this gives 



V < (-c+^K+lN{r,))5^\z{t)\) 



(22) 



+ [±K+%N{r 1 ) + l) 5{\z{t)\)\u{t)\ 
Setting x(s) = f ^(s/ 2 ). 11 follows from (21)-(22) that 
|u|oo < X (\x(t)\) |u|oo < x(2|*(t)|) 



Setting U H (a;, t) := V (x + R a (x, t), t), we see the deriva- 
tive V = V t (z, t) + Vt(z, t)z ofV(z, t) along (12) satisfies 

V = V t la] (x, t) + V™ ( x , t) {f(x, t, at) + u(i)}. (24) 

We deduce from (5), (21), and (23) that when the con- 
stant a (and so also 77) is sufficiently large, 

\u\< X (\x\) => 

Vt\x,t) + V^(x,t) [f(x,t,at)+u] < -§5 2 (\x\/2) 

and Si(\x\/2) < V™(x,t) < 5 2 (\x\ + 2S(\x\/2)/ V ) by 
(20). It follows that V la] is an ISS Lyapunov function for 
(12), so (12) is ISS for large a, by Lemma 6, as claimed. 
The UGAS conclusion is the special case where u = 0. To 
prove the iISS assertion, we instead follow the preceding 
argument up through (22) (which is valid since that part 
did use the unboundedness of S) and then substitute 
<f(|*(t)|)|u(t)| < cS 2 (\z(t)\)/2 + \u(t)\ 2 /(2c) into (22) 
and bound the resulting coefficient of \u(t) | 2 to show that 
V [a] is an iISS Lyapunov function for (12) for sufficiently 
large a (by again using (20) and (21) and taking v(s) = 
S 2 (s/2)c/4 and A(s) = fs 2 for a suitable f > 0), which 
implies that (12) is iISS for large a, by Lemma 6. 

We turn next to the special case where (2) is UGES. Let 
V satisfy the requirements of Lemma 1 above for (2), 
and let x(t) be any trajectory for (12) for any control 
u e U starting at x(t ) = x a . Define z(t) by (16). Ar- 
guing as before except with this new choice of V shows 
\R a (x{t),t)\ <2K\x(t)\/r) (by taking S(s) =Ks in (20)) 
and that (24) satisfies (by P3 and (13)) 

F<-|^)| 2 + c 3 |^)||u(i)| 

+c 3 \z(t)\ (IK 2 + %N(vj) {\x(t)\ + \u(t)\} 

< (-l+lc 3 K 2 + B V c 3 N(r])j \z(t)\ 2 

+ C3 |z(t)|(|^ 2 + fiV(r / ) + l)|u(t)|, 

since \x(t)\ < 2\z(t)\ for large 77 as before. If we now de- 
fine x € /Coo by x(s) = S /{8(1 + C3)}, then we deduce as 
in the UGAS case that if 77 is large enough, and if |u|oo < 
X(|a;(f)|) for all t, then we also have |u|oo < x(2|z(t)|) 
for all t and V < -~\z(t)\ 2 /2 < ~V{z(t),t)/{2c 2 ). This 
gives V(z{t),t) < V(z{t ),t )e 



-(t-t )/(2c 2 ) 



so 



ci 



^\x(t)\ 2 < ci\z(t)\ 2 < V(z(t),t) < c 2 \z(to)\ 2 e- — 
so our estimate on \R a (x(t), t)\ and the form of z(t) give 
\x(t)\ < ^(i + ^^le-^. (25) 



V < 



— + -K 

4 r/ 



■ V N( V )) 5 2 (\z(t)\)- 



We conclude as before that if (2) is UGES, then, when 
the constant a > is large enough, (1) is also UGES 



■5 



and (e.g., by the proof of l|Sontag fc Wandfl995l Lemma 
2.14)) (12) is ISES, which proves our theorem. 

Remark 10 The method we used in the proof of Theo- 
rem 7 can be used to prove the ISS property for (3) under 
appropriate growth assumptions on the matrix-valued 
function j:M"xIxR^r xm . Clearly, some growth 
condition on g is needed and linear growth of g is not 
enough, since x = —x + xu is not ISS. One way to extend 
our theorem to (3) is to add the hypotheses that g is C 1 
and that there is a constant c G > 1 such that for all t G M, 
x £ M™, and a > 0, \\g(x,t,at)\\ < c a + (6(\x\/2)) 1 / 2 , 
where 1 1 • 1 1 is the 2-norm on R nxm and 5 £ /C^ satisfies 
P1-P3 for some Lyapunov function V for (2). Applying 
the first part of the proof of Theorem 7 except with the 
new /Cqo function 



X(s) = 



cS(s/2) 



±{c + ^/5(Jj2)y 



(26) 



we then conclude as before that (3) is ISS for sufficiently 
large a > 0. If instead 5 £ JC is bounded, then (3) is iISS 
when a is sufficiently large, by our earlier argument. 

Remark 11 The decay requirement (5) on N £ M. 
from Theorem 7 can be relaxed, as follows. We assume 
the flow map of (2) satisfies the UGES conditions (6)- 
(7) for some D > 1 and A £ (0,K), where K satisfies 
(14), and we let V be as in Lemma 1. We can choose the 
constant C3 in (8) to be 

C3 = 4 fi-~A) 1) ' where Q = ^ K/X ~ 1 > ( 2? ) 

by the proof of Lemma 1 . It follows from our argument 
above that Theorem 7 remains true for cases where (2) 
is UGES if (5) is relaxed to 

¥>0s.t. sup V N( V ) < = • (28) 

n>n* 1LD(8 - l)c 

A similar relaxation can be made in the more general 
UGAS setting covered by Theorem 7. 

5 Proof of Theorem 9 

We first set U [a] = U [ t a] {x, t) + U™(x, t)f(x, t, at) for all 
x £ W l , t > 0, and a > 0. Using Assumptions HI. and 
H3. and (17) with p(t, I) independent of t gives 

f/H < -W(x,t)+p(at)G(x,t)-p(at)e(x,t) 

-(lL{r s P^)dl)d S )(^f+^) (29) 



< 



W(x,t)+ f*_ lP (al)dl \W{x,t) 
III (r s P( al ) dl ) ds\\W{x,t) 



along trajectories of (1), where we omit the argument 
(x, t, at) of /. For any a > 0, t > 0, and s £ [at — a, at] , 
Assumption H2. gives 



p(l)dl 



■(s) 



p(l)dl 



r(at) 



p(l)dl, 



(30) 



where t(u) := min{fcT : k € "L,kT > u} and 
t(u) := max{fcT : k G Z, kT < u}. The proof of 
(30) uses the facts that both endpoints of the interval 
[min{f (s), z(at)}, max{f (s), r(ai)}] are integer multi- 
ples of T and that, according to Assumption H2, the 
integral of p over any interval whose endpoints are inte- 
ger multiples of T is equal to zero. Choosing p max to be 
any global bound on \p(l) \ over all I £ K, (30) gives 



p(l)dl 



< 2Tp max y.se[at-a,at] (31) 



for all i > and a > 0. Hence, for all t > 0, 



p(al)dl 



t-i 



p(al)dl I ds 



1 

Q 
1 

cv 2 



p(l)dl 



at— a 
at 



at— a 



2T Pn 



1 

< — sup 

sd\at— a, at] 

< 2Tp max /a 



< 

p(l),n j ,/.s 

at 



P (l)dl 



(32) 



Using these estimates, (29) gives U lo,] < —W(x,t)/2 for 
all x £ K" and t > 0, as long as a > 8Tp max /c. Since 
W £ UPD, this gives the Lyapunov decay estimate. By 
Assumption H3. and (32), U [a] £ UPPD for large enough 
constants a > 0. This proves the theorem. 

6 Illustrations of Theorem 7 

We ne xt illustrate h ow Theorem 7 extends the re - 
sul ts of iKhalil (120021) and iPeuteman k AeveliU l|2002ft . 
In IPeuteman fc Aevelsl (|2002^ . the limiting dynamics 
(2) are assumed to be uniformly locally exponentially 
stable. In our first example, (2) is UGAS but not nec- 
essarily UGES. We next consider a class of systems 
(1) from IPeuteman fc Aevelsl 1)200 2') from identification 
where the limiting dynamics (2) is linear and exponen- 
tially stable. For these syst ems, our work complements 
IPeuteman fc Aevelsl (J2Q02) by providing formulas for 
Lyapunov functions for (1) that are expressed in terms 
of the quadratic Lyapunov functions for the limiting dy- 
namics and that have the additional desirable property 
that they are also ISS Lyapunov functions for (3) for 
suitable functions g. Finally, we apply our results to a 
friction model for a mass-spring dynamics. In all three 
examples, the limiting dynamics has a simple Lyapunov 
function structure so our results give explicit Lyapunov 
functions for the original rapidly time- varying dynamics. 
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6.1 Application to a UGAS dynamics that is not UGES 

Consider the following vari ant of the scalar example on 
l)Peuteman fc AevelsLl200l p.53): 



x = fix, t, at) = 

-cr x (a:)[2 + sfr±(f + cos(cr 2 (a;)))]{l + lOsin(crf)} 



(33) 



where a\ , a 2 : M — > M are C 1 functions such that a\ is 
odd, sup{|(ri(a;)| + \o-\{x)o-' 2 (x)\ : x £ M} < oo, <ri S /C 
on [0, oo), and 0" (s) < for all s > 0. One easily verifies 
the hypotheses of Theorem 7 using 

f(x,t) := -<Ji(x)[2 + sin(t + cos(a 2 (a;)))], 

V(x,t) = V(x) := / *oi(*)d*, 

<5(s) := 33<7i(2s), and iV(jj) := 60/ry 2 for large 77. 

This allows e.g. 01 (s) = 02(3) = arctan(s) in which 
case (2) is UGAS but not UGES because \x(t)\ < 2ir 
along all of its trajectories x(t). Condition Pi follows be- 
cause 01 (2s) < 2cti(s) for all s > 0, which holds because 
c" (s) < for all s > 0. Theorem 8 then gives the fol- 
lowing iISS Lyapunov function for (12) for large a > 0: 



i (0 /* sin(aZ) d/ ds 



(34) 



where := 2 + sin(Z + cos(ct2(C)))- in particu- 

lar, this is a Lyapunov function for x = f(x,t,at), 
and it is also an ISS Lyapunov function for (12) if 
S 6 K-oo (e.g. if (Ji (s) = sgn(s) ln(l + |s|) for |s| > 1 and 
02 (s) = arctan(s)). Our conditions on the 0i's cannot 
be o mitted even if the limiting dynamics (2) is UGES; 
see ijPeuteman fc Aevelst l2002t §8.2). For example, if 
cti(x) = x and a 2 {x) = x 2 , then (2) is UGES, but (33) 
is only sho wn to be locally expon e ntially stable for large 
a > 0; see iPeuteman fc Aevelsl l|2002() . This does not 
contradict our theorem because in that case (14) would 
be violated. 

6.2 A system arising in identification 

Consider the following v ariant of the example in 
l|Peuteman fc Aevelsl [20021 Section 8.1): 

x = f(at)m(t)m T (t)x + g(x,t,at)u, (35) 

with state x € M™ and inputs u € ffi m , where we assume 

i. / : M — > M is bounded and continuous and admits a 
o(s) function M and a constant /* < for which 

(1) /* = lim T ^ +00 27 J T T f(s)ds 

(2) I J t * 2 [f(s) - f*]ds\ < M(h - h) if t 2 > h 



ii. m : M. — > K™ is continuous and satisfies |m(t)| = 1 
for all i € R, and there exist constants a',/3',c > 
such that for alH € M, oc > 0, and x £ K™, we have: 

(1) oil < J t t+ " m(T)m T (T)dT < P'l 

(2) \\g(x,t,at)\\</3'{l + ^\x]} 

iii. j : 1" x 1 x R -> R" xm is continuous and is C 1 
in x, and there exists a constant K > 1 such that 
(d^Or, t, at)/dx\ < K Vx £ E™, f > 0, and a > 0, 
and each component of g. 

Here / denotes the n x n identity matrix. Notice that 
we allow / to take both positive and negative val- 
ues. T he special case of (3 5 ) whe re g = was stud- 
ied in IPeuteman fc Aevelsl l|2002f l where it is shown 
that the corresponding rapidly varying dynamics 
x = f(at)m(t)m T (t)x satisfies the hypotheses of Theo- 
rem 7 with the UGES dynamics 



x = f(x, t) := f*m{t) m T (t) x 



(36) 



and with S of the form 5(s) — rs for a constant f > 0. 
The particular case of (35) in which x = —m(t) m T (t) x 
has been extensively studied i n systems identification; 
see IPeuteman fc Aevelsl (|2002fl . However, these earlier 
results do not provide explicit ISS Lyapunov functions 
for (35). On the other hand, the following lemma pro- 
vides an explicit Lyapunov function for (36): 

Lemma 12 Let assumptions i.-iii. hold. If we choose 
P(t) = Kl+ [ [ m(l)m T (l)dlds, (37) 

Jt-c J s 

where n = c/(2\f*\) + -±rc 4 \f*\ thenV(x,t) =x T P(t)x 
is a Lyapunov function for ( 36) for which 2V/ a' satisfies 
the requirements of Lemma 1. 

To prove Lemma 12, we apply (17)-(18) with r = c and 
pit, I) = m(/)m T (/) and group terms to check that the 
derivative of V along trajectories of (36) satisfies 

V = i2f*K + c)x T mit)m T {t)x 

+2f*x T \St_ s J*m{l)m T (l)6lda\ mit)m T it)x 
—x T f.~ mil)m T il) dl x 

V < i2f*K + d)\m T it)x\ 2 -a'\x\ 2 

+2\r\\x\ [Sl_.jl\miirms] \mit)\\m T {t)x\ 
< (2/*/s + c)\m T it)x\ 2 + c 2 \f*\ \x\\m T it)x\ - a'\x\ 2 

everywhere, by ii-(l). By the triangle inequality, 



<?\f*\\x\\m T (t)x\ < ^a'lxl 2 + J-cV*\ 2 \™ T it)x\ 2 
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so V < (2f*K + c+^c 4 \f*\ 2 )\m T (t)x\ 2 ~ ±a'\x\ 2 
holds everywhere. Recalling that \m(t)\ = 1 everywhere 
and that P is everywhere positive definite, the fact that 
2V/a' satisfies the requirements of Lemma 1 follows 
immediately from our choice of k. Remark 10 now gives: 

Corollary 13 Let (35) satisfy conditions i.-iii. above 
and let V be as in Lemma 12. Then there exists a con- 
stant a > such that for each constant a > a , 



V [a] (x,t) :=V 



I - 



D(l)dl ds 



x, t 



where D(l) = (f (al) — f*)m(l)m T (I) is an ISS Lyapunov 
function for (35). 

6.3 Friction example 

The following one degree-of-f reedom mass-spring system 
from lde Queiroz et all l)2000j) arises in the control of me- 
chanical systems in the presence of friction: 

±i = x 2 

x 2 = — a i (at) x 2 — k(t)xi + u (38) 
~{a 2 (cd) + a 3 (at)e~^^} s&t(x 2 ) 

where x\ and x 2 are the mass position and velocity, re- 
spectively; <Xj, i — 1,2,3 denote positive time- varying 
viscous, Coulomb, and static friction-related coefficients, 
respectively; (3\ is a positive constant corresponding to 
the Stribeck effect; //(•) is a positive definite function also 
related to the Stribeck effect; k denotes a positive time- 
varying spring stiffness-related coefficient; and sat(-) de- 
notes any continuous function having these properties: 

(a) sat(0) = 0, (b) f sat(£) > V£ € M, 

(c) lim sat(0 = +1, (d) lim sat(f) = -1 ( 39 ) 

£ — >+oo £ — > — oo 

We model the saturation as the differentiable function 



sat(x2) = taxih(/3 2 x 2 ) , 



(40) 



where (3 2 is a large positive constant. Note for later use 
that |sat(x2) | < f3 2 \x 2 \ for all x 2 £ R. We assume the 
friction coefficients vary in time faster than the spring 
stiffness coefficient so we restrict to cases where a > 1. 

Our precise mathematical assumptions on (38) are: k 
and the o~i 's are bounded C 1 functions; [i has a globally 
bounded derivative; and there exist constants <7j, with 
fJi > and Gi > for i = 2,3, and a o(s) function 
s i ► M(s) such that 



(cTi(t) - <7i) di 



for all ti, t 2 £ R satisfying t 2 > t±. Although the cr,'s are 
positive for physical reasons, we will not require their 
positivity in the sequel. We show (38) satisfies the re- 
quirements of the version of Theorem 7 from Remark 10 
(with 6(s) = fs for a constant r) when (2) is 



X\ = x 2 
x 2 = -o\x 2 



[b 2 + a 3 e-^ x ^} s&t(x 2 ) - k(t) Xl 



(42) 



assuming this additional condition whose physical inter- 
pretation is that the spring stiffness is nonincreasing: 

3k Q , k>0 s.t. k a < k(t) < k and k'(t) < Mt > 0. 

To this end, set S := <T\ + (a 2 + cr 3 )(3 2 and 

V(x,t)=A(k(t)x\+xl)+x 1 x 2 , (43) 

where A := 1 + l/k a + [1 + S 2 /k ]/ax. Since Ak > 1, 
\(x 2 + x 2 ) < V(x,t) < A 2 k(\xi\ + \x 2 \) 2 for all x £ M 2 
and t > 0. Also, since k' < everywhere, the derivative 

V = V t (x,t) + V x (x,t)f(x,t) along trajectories of (42) 
satisfies 

V < V x (x, t)f(x, t) = [2Ak(t)xi + x 2 ]x 2 - [2Ax 2 + xi] 

x {a x x 2 + [d 2 + CTge-ft^ 2 )] sat(ar 2 ) + k(t)x x } 

and therefore, by grouping terms, we also have 

V < -k xj - (2Adi - \)x\ + S\x x x 2 \ (by (39)(b)) 

< ~b\x\ 2 - [!fx 2 + {Mi - l/2)x 2 - S\ Xl x 2 \] 

= -m 2 - k (Nil - £l^l) 2 + (£ + \ - ^i) 4 

< -b\x\ 2 , where b := min{fc D /2, Aai - 1/2}. 



< 



M(t 2 -h), i = 1,2,3 (41) 



The preceding inequalities imply that V/b is a Lyapunov 
function for (42) satisfying the requirements of Lemma 1. 
The integral bound requirement (13) from our theorem 
follows from (41) and the sublinear growth of tanh, since 
the integral bound can be verified term by term. The 
remaining bounds from (14) follow because \i and sat 
have globally bounded derivatives. We conclude that for 
sufficiently large constants a > 0, (38) admits the ISS 
Lyapunov function 

(£, t) = V Ui, 6 + ^ Jt^ /,* T a (l, dl ds, t\ 

where V is the Lyapunov function (43) for (42) and 

r a fl,0 ■■= {vi(al) + n a (l,£) tanh((3 2 b) (44) 

/i a (/,0 ~a 2 (al)- d 2 + (a 3 (al)-a 3 )e-^^ (45) 

so (38) is ISS for large enough a > 0, by Remark 10. 
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Remark 14 The preceding construction simplifies con- 
siderably if 02 and 03 in (38) are positive constants. In 
that case, the limiting dynamics (2) can be taken to be 

X\ = x 2 

x 2 = -oris*- {a 2 + aae-^^l sa,t(x 2 ) - k{t)x! 
and the ISS Lyapunov function for (38) becomes 

V (V, 6 + & JlWiial) - a x } dl d«) , tj 

with V defined by (43), since the /z Q (Z, £) tanh(/?2^2) 
term in the difference / — / is no longer present. 

7 Illustrations of Theorem 9 

Simple calculations show that Theorem 9 applies to (33) 
and (35) without controls (with the same choices of V 
that we used in our earlier discussions of those systems), 
assuming in the latter case that requirements i.-iii. hold 
and for instance / is a suitable periodic function and 
m(t) is bounded. We next show how Theorem 9 also 
applies to cases that are not tractable by Theorem 7. 

7.1 Dynamics that are not globally Lipschitz 

Simple calculations that we omit because of space con- 
straints show that the one-dimensional dynamics 

x = -x 3 + lOcos(crf) X „ (46) 
1 + ar 

satisfy H. using V(x, t) = x 4 /4, W(x, t)=x 6 , Q(x, t) = 
x 6 /(l + x 2 ), p(t) = 10cos(i), T = 2ir, and small enough 
c > 0, so (46) has the global Lyapunov function 

U [a] {x, t) = ^ - Qf (J 10 cos(al)c0j ds^j 

when a > is large enough. However, (46) is not covered 
by Theorem 7 since it is not globally Lipschitz in x. 

7.2 Dynamics with unknown functional parameters 
Consider the nonautonomous scalar system 

x 2 

x =p{at)—— + u , (47) 
1 + x z 

where p is an unknown, fast time- varying parameter sat- 
isfying \p(l)\ < a m for all I and some constant a m > 
and admitting a constant T > such that Assumption 



H2. holds for all k € Z. We assume the control u is ampli- 
tude limited in the sense that |tt| < u m for some constant 
u m > 0. We show that the saturated state feedback 

u = —u m arctan(x) (48) 

renders (47) UGAS. (A similar argument shows that u = 
—u m arctan(i?x) stabilizes (47) for any constant R > 1.) 
For simplicity, let a m — 10 and u m = 2. The derivative 
of V(x) = x 2 /2 along (47) in closed loop with (48) is 

V = -2x arctan(a;) + p(at) - 2 . (49) 

Simple calculations allow us to verify the hypotheses of 
Theorem 9 for the closed loop system using W(x,t) = 
2x arctan(x) , Q(x, t) = x 3 /(l + x 2 ), and small enough 
c > 0, so we know (48) indeed uniformly globally asymp- 
totically stabilizes (47) with control Lyapunov function 

U M { Xl t)= l -x 2 - (J* J' p(al) dl da) (50) 

when a > is sufficiently large. 

8 Conclusions 

Hypo theses similar to those of l|Peuteman fc Aevelsl 
2002, Theorem 3) are sufficient for uniform global 
exponential stability of rapidly time-varying nonlin- 
ear systems. We p ro vided complementary results to 
iPeuteman fc Aeve ls (2002) by establishing uniform 
global asymptotic stability of fast time-varying dy- 
namics without requiring local exponential stability 
of the limiting dynamics, and by constructing global 
Lyapunov functions for fast time-varying systems. Our 
Lyapunov constructions are new even in the special case 
of exponentially stable dynamics, and are input-to-state 
stable Lyapunov functions when the dynamics are con- 
trol affinc, under appropriate conditions. Our results 
apply to dynamics that are not necessarily uniformly 
Lipschitz in the state. We illustrated our methods using 
a friction example. 
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